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Abstract 

Recently Albert and Bousquet-Melou [2] obtained the solution to the long-standing 
problem of the number of permutations sortable by two stacks in parallel (tsip). Their 
solution was expressed in terms of functional equations. We show that the equally long¬ 
standing problem of the number of permutations sortable by a double-ended queue 
(deque) can be simply related to the solution of the same functional equations. Subject 
to plausible, but unproved, conditions, the radius of convergence of both generating 
functions is the same. Numerical work confirms this conjecture to 10 significant digits. 
Further numerical work suggests that the coefficients of the deque generating function 
behave as Kd ■ yC ■ where p = 8.281402207..., while the coefficients of the 

corresponding tsip generating function behave as Up ■ pA ■ n'*' with 7 « —2.473. The 
constants Kd and Kp are also estimated. 

Inter alia, we study the asymptotics of quarter-plane loops, starting and ending at 
the origin, with weight a given to north-west and east-south turns. The critical point 
varies continuously with a, while the corresponding exponent variation is found to be 
continuous and monotonic for a > — 1 / 2 , but discontinuous at a = — 1 / 2 . 


1 Introduction 

The problem of pattern-avoiding permutations appears to have been first considered by 
MacMahon [8]. However Knuth [7] was the first to consider a number of classic data struc¬ 
tures from the point of view of the permutations they could produce from the identity 
permutation (or, equinumerously, which permutations could produce the identity permu¬ 
tation). For the data structures considered - stacks and input-restricted deques - Knuth 
showed that of the n! possible input permutations of length n only Cn ~ const ■ 4” • 
and Sn rsj const - (3-|-2\/2)"’ • (re) could be sorted by stacks and input-restricted deques, 
respectively. This is a consequence of the fact that only 231-avoiding permutations are 
stack-sortable, as we discuss below. 
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Shortly thereafter, a number of authors, notably Evan and Itai [5], Pratt [9] and 
Tarjan [10] considered more general data structures. Foremost among these were two 
stacks in parallel, two stacks in series and deques. A deque, illustrated in Fig 2, is a double- 
ended queue, with insertions and deletions allowed at either end. Until recently [2], none 
of these had been solved. As remarked above, simple stacks can sort any permutation that 
does not contain three successive (but not necessarily consecutive) elements in the order 
231. We write this as the class Au(231), that is, the class of 231-avoiding permutations. 
For example 1573642 can’t be sorted by a stack as the elements 562 (among other sub¬ 
sequences) are in the forbidden relative order. Input-restricted deques are describable by 
the class Au(4231, 3241), whereas Pratt [9] showed that deques (without input or output 
restrictions) cannot be described similarly, as an infinite number of patterns would be 
needed in such a description. 

To establish a notation, let pn denote the number of permutations of length n that can 
be produced by two parallel stacks, let dn be the corresponding quantity for deques, and 
let Sn be the corresponding quantity for two stacks in series. We name the corresponding 
generating functions 

P{t) = D{t) = and S{t) = ^sX- 

In 2010, Albert, Atkinson and Linton [1] studied these problems with a view to estab¬ 
lishing upper- and lower-bounds to the relevant growth constants. For deques they found 
7.890 < Pd < 8.352 and for tsips they found 7.535 < ps < 8.3461, and commented that 
the actual growth constants may be equal, and may possibly be equal to exactly 8. As 
we show, the two growth constants do indeed appear to be equal, but to a slightly higher 
value, 8.28140... 

In 2015 Albert and Bousquet-Melou [2] found two coupled functional equations that 
give the generating function P{t). Unfortunately their representation does not allow for 
a single equation for the generating function, nor does it allow the asymptotics of pn to 
be obtained. However it does offer, in principle, a polynomial-time algorithm to obtain 
the coefficients pn- Other unanswered questions include the nature of the solution. Is it 
D-finite, or differentially algebraic? The answer to these questions is not known. 

In their solution, Albert and Bousquet-Melou first encoded the operations involved in 
sorting a permutation as words over the alphabet Ii, I 2 , Oi, O 2 , representing the input 
to, or output from, stack number 1 or stack number 2 respectively. Any sorting of a 
permutation can be effected by an operation sequence comprising a word in this alphabet, 
subject to certain constraints. They then pointed out that such words could be considered 
from two other points of view. The first is a mapping to quarter-plane random walks 
that return to the origin. More precisely, these are random walks in N x N that start and 
end at the origin. Making the identification Ii = N, I 2 = E, Oi = S, O 2 = W, where 
N, E, S, W, denote steps to the north, east, south and west respectively. Then operation 
sequences on words map to loops that return to the origin. The number of such loops of 
length 2n is given by CnCn+i, where Cn = Catalan number. However 

this is not a bijection as a given permutation may correspond to more than one loop. 

An alternative representation was given in terms of two-coloured arches, in which the 
operation sequences were encoded as arches drawn between numbered points on a line, 
with arches above the line being of a different colour than those below the line. For further 
details of these connections the reader is referred to [2]. 
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In this work our principal result is that the generating function for deques can be 
simply related to the generating function for two stacks in parallel. We therefore provide 
a comparable solution to the deque problem to that given in [2] for the problem of two 
stacks in parallel in Theorem 3.9, which states that the deque and tsip generating functions 
D and P, defined above, satisfy the following two (equivalent) equations; 

^ ’ 2t- {D{t) -1-t- D{t)) ^ ^ 

and 

2D{t) = 2 + t + 2tP{t) - 2t^P{t) - - AP{t) + 4P(t)2 - 8tP(t)2 + At'^P‘^{t) - 4tP(t). 

( 2 ) 

Subject to the veracity of conjectures 10, 11 and 12 in [2], we prove that the radii of 
convergence of the two generating functions P{t) and D{t) are the same. Further, with one 
extra assumption, in Theorem 4.6 we prove that D{t) = D{tc) + ko^/t — tc + o{y/t — tc) 
for some constant kjj 

For both these problems we also give a detailed (numerical) study of the asymptotics, 
based on a 500 term expansion of the generating functions. In this way we estimate the 
radius of convergence to 10 significant digits for both problems, and find agreement at 
that level of precision, thus strengthening our confidence in the aforementioned conjecture 
that the critical points are identical. We also find the leading and next-to-leading critical 
exponents, which, perhaps surprisingly, are different for the two problems. For two stacks 
in parallel we find 

P{t) = P{p) - t,P\p)il - t/tc) + Pl{tc){l - + P2(tc)(l - + o((l - 


where Pi(tc), P 2 {tc) etc. are constants that occur as the lowest order terms in the ex¬ 
pansion of implicitly defined functions Pi{t), P2{t) etc. around t = tc- We estimate 
tc ~ 0.1207524976, 71 « 1.473, and 72 ~ 1.946. Furthermore 

= 2{l-^c?- 

For deques we find 

D{t) = D{tc) + Pl(tc)(l - t/tcf/"^ + D2{tc){l - t/tcV^ + 0(1 - ^/^c)^^ 


where Di{tc), T> 2 (ic) etc. are constants that occur as the lowest order terms in the ex¬ 
pansion of implicitly defined functions Di{t), D 2 {t) etc. around t = tc-- We estimate 
tc ~ 0.1207524977, and 73 « 0.973. Furthermore 


D{tc) 


1 + t 


3/2 

C 


l-tc’ 


and 

D,{tc) = - 23/4 . ( 3 ) 

We have been unable to obtain a solution to the more difficult problem of the number 
of permutations of length n sortable by two stacks in series, as shown in Fig. 1, but have 
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Figure 1: Allowed input and output sequences for two stacks in series. 


obtained the most extensive enumerations to date of the coefficients Sn of the generating 
function S(t}, having obtained these exactly for n < 20 and approximately for longer 
permutations. An analysis of the asymptotics based on these enumerations is given in [4]. 

The associated problem of weighted quarter-plane loops also displays some unusual 
and hence interesting features for a simple lattice model. Recall that the loops can take 
steps in the ±x and ±y directions, and start and end at the origin (0,0), and that all 
vertices must have non-negative co-ordinates. There is, in addition, a weight, or fugacity 
a associated with NW and ES corners, so the generating function can be written 

Q{u, o) = ^ qn,mU^ ■ a™ 

where qn,m is the number of such loops of n steps with m weighted corners. We are 
interested in the asymptotic behaviour of the coefficients, which is expected to be 

[vP']Q{u, a) ~ const ■ qc{a)~'^ ■ 


As shown by Albert and Bousquet-Melou [2], not only does the critical point qdo) vary 
with a, which is to be expected, but so does the critical exponent g{a), which is unexpected. 
In particular, this does not happen for similarly weighted loops in the half-plane or full- 
plane. Based on our numerical studies, we conjectured the variation of qc{o) with o, 
providing numerical confirmation for the conjectured dependence in [2]. We also studied 
the critical exponents, and conjectured these for various values of a, but were not able 
to guess the a dependence. In private correspondence, Kilian Raschel has given us his 
conjectured value for the exponents which is 


9{a) 


arccos 


vr 


g-l \ ’ 


for a > 0. 


This agrees with our numerical results, not only for a > 0, but more broadly for a > —1/2. 
However for a = —1/2, we hnd the exponent takes the value 3/4, rather than 1 as given 
by the above formula. Thus we have the additional interesting feature of a discontinuous 
critical exponent arising in this simple lattice walk model. 

In the next section we establish our notation, define certain operations and prove a 
number of lemmas and propositions necessary for the proof of the main theorem. In the 
following section we derive the generating functions for various quantities, such as bi¬ 
coloured Dyck paths, weighted loops as defined above, and for a special class of sequences. 
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defined below. Putting these definitions together leads us to a functional equation con¬ 
necting D(t) and P{t), from which we obtain our principal result Theorem 3.9. In section 
4 we analyse this functional equation, and, subject to a conjecture in [2], prove Theorem 
4.6, mentioned above. In section 5 we first discuss the asymptotics of weighted quarter- 
plane loops, then that of deques and two stacks in parallel. It is only through a remarkable 
cancellation that these two generating functions have different critical exponents. For de¬ 
ques we conjecture a square-root singularity, while for tsips the exponent arising in the 
generating function is approximately 1.473, which does not suggest any obvious rational 
number with low denominator. In a wild speculation we give a possible exact value for 
this exponent. The final section gives our conclusions. 

2 Notation and canonical operation sequences 

Consider the operations Ii, l2,0i,02, where Ii and I 2 represent input to the top and 
bottom of the deque, respectively, and Oi and O 2 represent output from the top and 
bottom of the deque, respectively, as shown in Fig 2. We will call a permutation deque- 
achievable if it can be produced by a deque, and we will call it tsip-achievable if it can 
be produced by two stacks in parallel. Each achievable permutation can be encoded as a 
word over the alphabet {/i, 125 Oi, O 2 }, such that the total number of Is is the same as 
the total number of O’s, and any prefix contains at least as many I’s as Os. We call such 
a word an operation sequenee. For example, the permutation 4123 is achievable as it is 
produced by the operation sequence I 1 I 1 I 1 I 2 O 2 O 2 O 2 O 2 (see Fig 3). This permutation is 
not, however, tsip-achievable. 



Output 




Input 


Figure 2: The input and output operations Ii, I 2 , Oi and O 2 on a deque. 

We call an operation sequence w a tsip word if it has the additional properties that 
for j = 1,2 the total number of Ij’s is equal to the total number of Oj’s, and each prefix 
contains at least as many Ij’s as Oj’s. In this case, the operation sequence w corresponds 
to a tsip-achievable permutation. Note that if a permutation is tsip-achievable, then it is 
also deque-achievable, but the converse does not hold. 

Definition 2.1. We call an operation sequence w canonical if it has the following three 
properties: 
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Figure 3: The sorting procedure I 1 I 1 I 1 I 2 O 2 O 2 O 2 O 2 produces the permutation 4123. This 
is the shortest permutation which cannot be produced by two stacks in series. 

• (outputs eagerly) w contains no consecutive terms I 1 O 2 or I 2 O 1 . 

• (standard) Any tsip sub-word of w begins with Ii. 

• (top happy) Whenever one of the letters I 2 or O 2 appears in w, the number of 
preceding Is in w must be at least two greater than the number of preceding Os. 
In other words, these two letters only appear when there are at least two numbers 
already in the deque. 

Note that this definition also defines what is meant by the terms outputs eagerly, 
standard and top happy. 

In Proposition 2.5 we will prove that every deque achievable permutation vr is produced 
by a unique canonical operation sequence. This is analogous to the result in [2], which 
in our language states that every tsip achievable permutation vr is produced by a unique 
standard tsip-word which outputs eagerly. 

Lemma 2.2. Any achievable permutation vr can be produced by a canonical operation 
sequence. 

Proof. Assign the ordering Oi < O 2 < Ii < I 2 to the letters. Since tt is achievable, it is 
produced by some operation sequence. Let w be the operation sequence which produces 
TT and which is lexicographically minimal. Then we just need to show that w is canonical. 

If w contains either / 1 O 2 or / 2 O 1 as a sub-word, then we can switch these letters to 
produce w' < w, where w' also produces vr. But this contradicts the minimality of re, so 
w does not contain either I 1 O 2 or I 2 O 1 . Hence w outputs eagerly. 
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Suppose that w contains a tsip sub-word v which begins with 12- Then we can form 
v' by changing each subscript in v to the other subscript. Note that v' begins with Ii, 
so v' < V. Moreover, v and v' have the same effect on the state of the machine as each 
other. Hence, \iw = avb, then the word w' = av'b produces the same permutation vr. But 
w' < w, which contradicts the minimality of w. So w does not contain a tsip sub-word 
beginning with 12- Hence, w is standard. 

Now suppose that w can be split up as tc = uv, where the number of Is in u is at most 
one more than the number of Os, and where v begins with either I 2 or O 2 . Then we can 
form v' by changing each subscript in v to the other subscript. Note that v' begins with 
Ii if V begins with I 2 , and v' begins with Oi if v begins with O 2 , so in either case v' < v. 
After u is applied to the machine, the deque contains at most one element, so v' and v 
have the same effect on the machine, with v' essentially doing everything on the opposite 
end of the deque to v. Hence uv' produces vr. Moreover, uv' < uv = w, contradicting 
the minimality of w. Therefore, w cannot be split as uv in this way, so whenever one of 
the letters I 2 or O 2 appears in w, the number of preceding Is in w must be at least two 
greater than the number of preceding Os. Hence, w is top happy. 

Therefore, w is standard, top happy and outputs eagerly. Hence w is canonical. □ 

As in [2], we will consider the type of an operation sequence w: 

Definition 2.3. The type of an operation sequence w is the word over the alphabet {/, O} 
obtained by deleting all of the subscripts of w. 

Lemma 2.4. If two operation sequences u and v produee the same permutation vr and 
both output eagerly, then u and v have the same type. Moreover, the ith operation of u 
moves the same item as the ith operation of v. 

Proof. Let u = U1U2 • .. U2n and let v = V1V2 ... V2n- We will prove the lemma by induction. 
Assume that for some k, the words ui.. .u^ and vi.. .v^ have the same type. We will 
prove that the letters Uk+i and Vk+i also have the same type, and move the same item 
as each other. After the fcth operation, both u and v have the same number of input 
steps and the same number of output steps as each other, so, since they produce the same 
sub-permutation vr, the items currently in the deque according to u must be the same as 
the items currently in the deque according to v, though the order may be different. Since 
u and V output eagerly, they will both output at this point if and only if the next item to 
be output is currently in the deque. Since this is the same for both u and v, the letter u^+i 
and Vk+i have the same type. Moreover, if they are both output steps, then the element 
for each will be the next element to be output, whereas if they are both input steps, the 
item moved will be the next element in the input. In both cases the item moved is the 
same for both u and v. This completes the induction, so u and v have the same type. □ 

Proposition 2.5. Every deque-aehievable permutation vr is produced by a unique canonical 
operation sequence. 

Proof. We have already shown that every deque-achievable permutation is produced by 
some canonical operation sequence, so it remains to show that this operation sequence is 
unique. Let u and v be two canonical operation sequences that produce the permutation 
IT, with u < V. Then we just need to show that u = v. 
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Suppose for the sake of contradiction that u < v. Let u = ui...U 2 n and let v = 
• • • V 2 n, and let k be minimal such that Uk ^ v^, so < v^- Note then that before the 
kth move, the position of the deque, input and output is the same according to either u 
or V. By Lemma 2.4, u and v have the same type, so the letters and have the same 
type. Therefore, Uk,Vk are either Oi, O 2 , respectively, or Ii, I 2 , respectively. Since Vk is 
O 2 or I 2 , and v is top happy, there must be at least two items in the deque before this 
move. Now, by Lemma 2.4, the element moved by Uk and Vk must be the same, but Uk and 
Vk act at opposite ends of the deque, so they must both be input operations. Therefore, 
Uk = Ii and Vk = l 2 - We will now obtain a contradiction by showing that the operation 
sequence v is not standard. 

Define the sequences fco, ki,... ,km and Iq, li,... ,1m, of indices in {k, ..., 2n} as follows: 

• ko = k 

• For each i, the operations Vk^ and Vl^ move the same item, with Vk^ an input and vi. 
an output move. In particular, this determines Iq. 

• For i > 0, the values ki+i and k+i are chosen inductively so that ki < fcj+i < k < 
li^i. Moreover, they are chosen to maximise k+i. If no such /cj+i and Zj+i exist, 
then we set m = i and terminate. Note that the process must terminate at some 
point since li,l 2 ,... is a sequence of integers which is increasing and bounded above. 

For 0 < i < m, let ti £ {1, 2} be the value such that = If. 

Let D be the set of items which are in the deque just before operation k. Let Vj be the 
first operation which outputs an item in D, and let this item be d. Then d must be at one 
of the ends of the deque just before operation k. We will prove the following statements 
for 0 < i < m in a single induction: 

• k < ki < j. 

• Uki = h-u- That is, Uk^ and Vki put the element they move at opposite ends of the 
deque. 

• k < j- 

• u/. = Oti. So Vki and vi- act on the same side of the deque. 

• Vl^ = Oti. 

• ui^ = Os-ti. So Uki and Ul^ act on the same side of the deque. 

• ui. = Os-ti- 

For the base case, by the definition of k and j, we have ko = k < j and UkQ = Uk = 
h = h-tQ- 

Now for the inductive step. We will assume that k < ki < j and Uk^ = H-u and 
show that k < j, Vl^ = Ot^, ui^ = k < h+i < j and = h-u+i- Let a* be 

the item moved by operations , u;. and Vl^. Since . = h-ti and = Iti , we 

have {ufe-jUfe.} = {Ii,l 2 }. Therefore, one of Vk^ and Uk^ places a* on the end of the deque 
which d has to be removed from. So a* must be removed before move j, hence li < j. 
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Since there is an item d on the deque throughout the whole time that a* is on the deque, 
the input and output operations which move a* must act on the same side of the deque. 


Hence, u;. = Oj. and u/. = Oa-ti- Now to finish the induction, we just need to show that 
k < /cj+i < j and Ufci+i = h-ti+i, in the case where i < m. First, k < ki < /cj+i < k < j. 

Also, by construction, the output move vi^ = O*- for Oj occurs between the input 
and output moves for Ui+i, so we must have = Iz-u-, so tj+i = 3 — tj. Similarly, 

= hi = ■ This completes the induction. 

We will now show that w = VkVk+i ■ ■ -vi^ is a tsip word. Since Im < j, nothing is 
output from the deque from D. Now, suppose some item is input but not output by w. 
Let the input and output moves of that item be Va and Vb, so k < a < Im < b and let 
i be maximal such that ki < a. If i = m then km < a < Im < b, which means that we 
shouldn’t have terminated at m. If i < m then ki < a < fej+i < k < k+i < b, which 
contradicts the maximality of h+i- Hence every item which is input by w is also output by 
w. Therefore, Vk ■ ■ ■ vi^ is a tsip sub-word of v which begins with I 2 , so v is not standard. 
This contradicts the statement that v is canonical, so we must have u = v. Therefore, any 
two canonical operation sequences which produce the same permutation are equal. □ 

3 Enumeration 

In the previous section we showed that every deque-achievable permutation is produced by 
a unique canonical operation sequence, hence the number of deque achievable permutations 
of length n is equal to the number of canonical operation sequences of length 2n. Recall 
that this sequence is given by the generating function D{t). 

Before we proceed any further, we will define the generating functions Q and P in the 
same way as in [2]. Let Q{a, u) be the generating function for quarter plane loops, counted 
by half-length, conjugate to the variable u and number of NW or ES corners, conjugate 
to a. Let P{t) be the length generating function for permutations which can be produced 
by two stacks in parallel. We note that this is the same as the generating function for 
standard tsip words which output eagerly, counted by half-length. Now using corollary 9 
from [2], we have the following relation between these generating functions: 



(4) 


Now we will define some new generating functions to characterise the length generat¬ 
ing function D{t) of deque-achievable permutations. We will call a non-empty operation 
sequence unbreakable if it contains no tsip-sub-words, apart from possibly itself, and the 
deque is never empty during the associated procedure. Let M{a,u,x) be the generating 
function for top happy, unbreakable operation sequences, counted by half-length, conju¬ 
gate to u, number of appearances of a sub-word I 1 O 2 or I 2 O 1 , conjugate to a, and number 
of times when there is only one element in the deque, conjugate to x. 

We define a bi-coloured Dyck path to be a Dyck path where each step is coloured red or 
blue. A multicoloured peak is an up-step followed by a down-step of the opposite colour. 
Let T{a,u,x) be the generating function for bi-coloured Dyck paths, where every step 
from height 0 or I is red, counted by half-length, conjugate to n, number of multicoloured 
peaks, conjugate to a, and number of vertices of the path at height 1 , conjugate to x. 
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We define a new generating function Qi{a,u,x) for quarter plane loops where every 
step from the origin is a north step, counted by half-length, conjugate to u, number of 
NW or ES corners, conjugate to a, and number of steps to the origin, conjugate to x. This 
is related to Q by Q{a,u) = Qi{a,u,2). We can also think of Q and Qi as generating 
functions for tsip words, by replacing each N, E, S, W step in the quarter plane loop with 
/i, I 2 , Oi, O 2 , respectively. Then Q{a,u) is the generating function for tsip operation 
sequences, counted by half-length, conjugate to u and number of appearances of a sub¬ 
word I 1 O 2 or hOi, conjugate to a. The generating function Qi{a,u,x) is similar, except 
that it only counts those tsip words where every letter which appears after an equal 
number of Is and Os is an Ii. The other difference is that Qi also counts these words by 
the number of times there are an equal, non-zero number of Is and Os, conjugate to x. 


Lemma 3.1. The generating function for bi-coloured Dyck paths T{a, u, x) is given by the 
equation 


T{a, u, x) 


4 -|- 2xu — 2xau — x xVl — 12u -|- 4u‘^ — Aau — 8au‘^ 

4 — 2xu — 2xau — x x\/Y^^T2fri^4u?~^^^afffl^4ahE~^^8a^ 


(5) 


Proof. Let Ti{a,u,y) be the generating function for bi-coloured Dyck paths, where each 
step from height 0 is red, counted by half-length, conjugate to u, multicoloured peaks, 
conjugate to a, and number of vertices of the path at height 0, including the end points, 
conjugate to y. Let T 2 {a, u) be the generating function for bi-coloured Dyck paths, counted 
by half-length, conjugate to u, and multicoloured peaks, conjugate to a. Note that each 
multicoloured Dyck path cj satisfies exactly one of the following: 


• a; is empty, 


• a; is made up of a peak followed by a (possibly empty) multicoloured Dyck path. 


• a; is made up of an up-step, followed by a non-empty multicoloured Dyck path, 
followed by a down-step, followed by a (possibly empty) multicoloured Dyck path. 

The contribution to T 2 from the first case is simply 1. 

In the second case, the 4 possible peaks are counted by u{2a + 2), so the total contri¬ 
bution to T 2 from this case is u{2a + 2 )T 2 . 

In the third case, the contribution is simply 4ri(T2 — 1)T2. The multiplier 4 arises 
because there are four possibilities for the colours of the up step and down step. T 2 — 1 
counts the first (non-empty) path and T 2 counts the second path. The term u is in the 
expression because the half length of oj is one more than the sum of the half lengths of 
the shorter paths. Hence we get the equation 


T 2 {a, u) = 4rt(r2 - l)r 2 -|- u{2a -|- 2 )T 2 -|- 1. 


Solving the quadratic gives 

1 -|- 2tt — 2au ± Vl — 12u -|- 4u‘^ — 4au — 8au‘^ 


T2{a,u) = 


8u 


Note that since T 2 {a, 0) = 1, we must use the negative square root, so 

l + 2u — 2au — a/1 — 12n -|- — Aau + la^u^ — 8au^ 


T2{a,u) = 


8u 
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We now calculate Ti(a, u, y) in a similar way. Each multicoloured Dyck path u, where 
each step from height 0 is red, satisfies exactly one of the following; 

• a; is empty, 

• a; is made up of a peak followed by a (possibly empty) multicoloured Dyck path, 

• a; is made up of a red up-step, followed by a non-empty multicoloured Dyck path, 
followed by a down-step, followed by another (possibly empty) multicoloured Dyck 
path where each step from height 0 is red. 

The contribution to Ti from the first case is simply y. 

In the second case, the two possible peaks are counted by uay + uy, so the contribution 
to Ti from this case is {uay + uy)Ti. 

In the third case, the half length of oj is one more than the sum of the half lengths of the 
shorter paths, the number of vertices at height 0 is one more than the number of vertices 
at height 0 in the second path. Hence the contribution from this case is 2uy{T2 — l)Ti. 
Hence, we get the equation 

Ti (a, u, y) = 2uy{T2 - l)Ti {uay + uy)Ti + y. 


Solving this using the formula for T 2 gives 
Ti{a,u,y) = 


4y 


4 -|- 2yu — 2yau — y + yVl — 12u -|- — Aau + — Sau^ 

Finally, we can calculate r(a, u, x) in a similar way. Each multicoloured Dyck path w, 
where each step from height 0 or 1 is red, satisfies exactly one of the following; 


• a; is empty 

• cj is made up of a red up-step, followed by a (possibly empty) multicoloured Dyck 
path, where each step from height 0 is red, followed by a red down-step, followed by 
another (possibly empty) multicoloured Dyck path where each step from height 0 or 
1 is red. 


The contribution to T from the first case is simply 1. In the second case, the half length 
of (jO is one more than the sum of the half lengths of the shorter paths, and the number 
of multicoloured peaks in uj is equal to the sum of the numbers of multicoloured peaks in 
the two shorter paths. The number of vertices at height 1 in the long path is equal to the 
number of vertices at height 0 in the first short path, plus the number of vertices at height 
I in the second short path. Hence, 

r(a, u, x) = uTi{a, u, x) ■ T{a, u,x) + 1 

Solving this using the formula for Ti gives 

^, .4-1- 2xu — 2xau — x + xVl — 12u + — Aau + — Sau^ 

1 [a,u,x) = - ^ 

4 — 2xu — 2xau — x + xv 1 — 12u + Au'^ — Aau -|- 40 ^^^ — Sart^ 

as required. □ 
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Lemma 3.2. The generating function Qi{a,u,x) is given by the equation 


Qi{a,u,x) 


2Q 

2Q — xQ + X 


(6) 


Proof. Let U{a, u) be the generating function for non-empty quarter plane loops, which 
only touch the origin at the start and end, counted by half-length, conjugate to u and 
number of NW or ES corners, conjugate to a. Then each loop counted by Q can be 
written uniquely as a sequence of loops counted in f7, so we get the equation 


Q{a,u) 


1 

1 — U{a, u) 


Or, equivalently, 


U{a,u) = 1 — 


1 

Q{a,u)' 


Now, by reflecting a loop about the line x = y, we see that amongst the loops counted 
by U with a given half-length and number of NW or ES corners, exactly half begin with 
a north step. Hence, each loop counted by Qi can be written uniquely as a sequence of 
loops counted in . Moreover, the power of x in the corresponding monomial in Qi is 
equal to the number of terms in the sequence of loops from Therefore, we have the 
equation 


Qi{a,u,x) 


1 

1 — \U{a, u)x 


Hence, we can write Qi in terms of Q as desired 


Qi{a,u,x) 


1 

2Q 

2Q — Qx + X 


□ 


Before we can relate these generating functions to M and D, we will need to consider 
a decomposition of operation sequences, as described in the lemma below. In lemma 3.5 
we will show that this decomposition is unique. 

Lemma 3.3. Given a non-empty operation sequence w it is possible to decompose w as 


W = XIWIX2W2 . . . X2m-lW2m-lX2mV, 

where m G Z>o, each Xi G {Ii, Oi, / 2 ) O 2 }, such that X 1 X 2 ■ ■ ■ X 2 m is an unbreakable opera¬ 
tion sequence, each Wi is a (possibly empty) tsip sub-word, and v is an operation sequence. 

Proof. We construct the decomposition as follows: First we decompose w as w = uv, so 
that the first point at which the deque is empty is immediately after u is applied, in other 
words, u is the minimal, non-empty prefix of w which is also an operation sequence. Note 
that we may have u = w. Since the deque is empty before and after v is applied, v is an 
operation sequence. 
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Now, xi must be the first letter of u. Then, for each i, we let Wi be the longest tsip- 
sub-word of u which starts immediately after Xi. Note that this necessarily exists, since 
Wi can always be the empty word. Then Xj+i is the letter after Wi. We continue this 
until we have decomposed all of u as u = xiwi... Wj-iXjWj. Note that u = xiwi.. .Wj 
and Wj both contain an equal number of Is and Os, so xiwi.. .Xj also contains an equal 
number of Is and Os. But since u is minimal, Wj must be empty. Note that by definition, 
each Wi is a tsip-sub-word. So we just need to prove that xiX 2 ■ ■ - xj is an unbreakable 
operation sequence. Since xiwi... Wj-iXj is an operation sequence, and each Wi is an 
operation sequence, it follows that xiX 2 ■ ■ ■ xj is an operation sequence. Also, if xiX 2 ■ ■ - Xi 
is an operation sequence for any i < j, then xiwi .. .Xi is also an operation sequence, 
which contradicts the minimality of u. Now we just need to show that xi... Xj contains 
no tsip sub-words, other than itself and the empty word. Suppose that it does contain 
some tsip sub-word Xa ■ ■ - Xi, with i > a. If a > 1 then Wa-iXaWaXa+i ■ ■ - Xi is also a tsip 
sub-word, which contradicts the maximality of the length of Wa-i- If a = 1 and i < j, 
then XaXa+i ... x* is not an operation sequence, so it is certainly not a tsip sub-word. 
Hence a = 1 and i = j, so xiX 2 ■ ■ - Xj contains no tsip sub-words other than possibly itself. 
Therefore xiX 2 ■ ■ - Xj is unbreakable. □ 

Lemma 3.4. Let w be a non-empty operation sequence with the decomposition 


W = X 1 W 1 X 2 W 2 . . . X2m-lW2m-lX2mV, 

where m G Z>o, each x* G {Ii, Oi, I 2 , O 2 }, such that X 1 X 2 ... X 2 m is an unbreakable opera¬ 
tion sequence, each Wi is a (possibly empty) tsip sub-word, and v is an operation sequence. 

Then any tsip sub-word of w, which is not a prefix, is contained in one of the words Wi or 

V. 

Proof Let u = X 1 W 1 X 2 W 2 ■ ■ ■ X 2 m-i'W 2 m-iX 2 m- Suppose for the sake of contradiction that 
some tsip sub-word w' of w is not contained in any Wi or v, and w' is not a prefix of w. 

First we consider the case where w' does not intersect with v. Then let w' = sxiWiXi+iWi+i... xjt 
be a longer tsip sub-word of u, where f is a prefix of Wj and s is a suffix of Wi-i. Since 
f is a suffix of an operation sequence w' and a prefix of another operation sequence, it 
must be an operation sequence. Moreover, since t is a prefix of a tsip word, t is also a 
tsip word. Similarly, s is a tsip word. Hence, s,Wi,Wi^i, ... ,Wj-i,t are all tsip words, 
as is w' = swiXi+iWi+i ... Wj-iXjt, so Xi+iXi +2 ■ ■ - xj is also a tsip word. But this is a 
contradiction, since X 1 X 2 ... X 2 m is unbreakable. 

Now we consider the case where w' intersects with v. Then we can write w' = 
where u' is a suffix of u and v' is a prefix of v. Since w' is not a prefix of w, the word u' 
is also not a prefix. Since w' is not contained in v, the word u' is non-empty. Since u' is a 
prefix of an operation sequence and a suffix of an operation sequence, it is also an operation 
sequence. But then u' satisfies the properties of w' in the first case, a contradiction. □ 

Lemma 3.5. Given a non-empty operation sequence w there is a unique way to decompose 
w as 

W = X 1 W 1 X 2 W 2 ■ . . X2m-lW2m-lX2mV, 

where m G Z>o, each x* G {Ii, Oi, I 2 , O 2 }, such that X 1 X 2 ... X 2 m is an unbreakable opera¬ 
tion sequence, each Wi is a (possibly empty) tsip sub-word, and v is an operation sequence. 
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Proof. Let w = X 1 W 1 X 2 W 2 ■ ■ ■ X 2 m-iW 2 m-iX 2 mV be one such decomposition of w. We just 
need to show that this is the same as the decomposition constructed in lemma 3.3. 

First we will show that u = X 1 W 1 X 2 W 2 ■ ■ ■ X 2 m-iW 2 m-iX 2 m is the shortest operation 
sub-sequence starting from the start of w. Since each Wi is an operation sequence, 
and xi... X 2 m is an operation sequence, X 1 W 1 X 2 W 2 ■ ■ ■ X 2 m-iW 2 m-iX 2 m is an operation se¬ 
quence. Let u' be a non-empty prefix of u such that u' / u. Then u' = xiwi... Xi-iWi-iXit 
for some i < 2m, where t is a prefix of Wi. Then, since xi... X 2 m is unbreakable, there are 
strictly more Is than Os in xi... x,. Also, each Wk contains an equal number of Is and 
Os and t contains at least as many Is as Os. Hence u' contains more Is than Os, so it is 
not an operation sequence. Therefore, u is the shortest non-empty prefix of w, which is 
also an operation sequence. 

Finally, by the previous lemma, every tsip sub-word of w is either a prefix of w or 
is contained in one of the words Wi or v. Hence, each wt is the tsip sub-word of max¬ 
imal length with that starting point. Therefore, this decomposition is the same as the 
decomposition constructed in lemma 3.3 □ 

Lemma 3.6. Letw be a non-empty operation sequence, with decomposition w = xiwi... W 2 m-iX 2 mV. 
Then w is canonical if and only if the following conditions hold: 

• The (unbreakable) operation sequence s = X 1 X 2 ■ ■ ■ X 2 m is top happy, 

• Each tsip word Wi is standard and outputs eagerly, 

• V is canonical, 

• If some XiXi+i is either I 1 O 2 or I 2 O 1 , then Wi is non-empty. 

Proof. If w is canonical, then w is top happy, so s is top happy and v is top happy. Since 
w is standard and outputs eagerly, any sub operation sequence of w is also standard and 
outputs eagerly; in particular this includes each Wi, as well as v. Hence v is canonical. 

Finally the fourth condition follows immediately from the fact that w outputs eagerly. 

Now we will assume the four conditions and prove that w is canonical. A sub-word of 
the form 10 can only appear in w inside one of the sub-words Wi or v or as XjXj+i, where 
Wi is empty. The conditions clearly make it impossible for such a sub-word to be either 
I 1 O 2 or I 2 O 1 , so w outputs eagerly. Now we consider re to be a bi-coloured Dyck path. 

Since s and v are top happy, every step from height 0 or 1 in u) which comes from s or u 
is red. Any step from height 0 or 1 which comes from some Wi, must be at height 1 in re 
and height 0 in Wi. Then this step must be red since Wi is standard. Finally, by Lemma 
3.4, any tsip sub-word of w is either a prefix of w, contained in v or contained in one of 
the words Wi. Hence, since each Wi and v are standard, and w begins with Ji, the word w 
is also standard. Therefore, w is canonical. □ 

Now recall that the generating function M (a, u, x) is the generating function for top 
happy, unbreakable operation sequences, counted by half-length, conjugate to u, number 
of appearances of a sub-word I 1 O 2 or I 2 O 1 , conjugate to a, and number of times when 
there is only one element in the deque, conjugate to x. 


14 


Lemma 3 . 7 . The generating function M satisfies the equation 

T{a, u,x) = M (^1 + ^ ^ 

Proof. We first note that given an operation sequence w, there is a corresponding bi¬ 
coloured Dyck path, formed by replacing each Ji with a red up-step, each I 2 with a blue 
up-step, each Oi with a red down-step, and each O 2 with a blue down-step. Then the 
condition that the steps from height 0 or 1 in the Dyck paths counted by T are all red 
is equivalent to the condition that the corresponding operation sequence is top happy. 
Hence, we can consider T{a,u,x) to be the generating function for top happy operation 
sequences, counted by half-length, conjugate to u, number of consecutive steps I 1 O 2 or 
I 2 O 1 , conjugate to a, and number of times in the procedure that the deque contains exactly 
one element, conjugate to x. 

For each non-empty operation sequence w counted by T{a,u,x) we consider the de¬ 
composition w = X 1 W 1 X 2 W 2 ■ ■ ■ X 2 m-iW 2 m-iX 2 mV described in lemma 3.5. In particular, 
we consider the contribution to T from all words w with a given unbreakable operation 
sequence s = X 1 X 2 ■ ■ ■ X 2 m- Since w is top happy, s is also top happy, so we will only con¬ 
sider the top happy, unbreakable operation sequences s, which are exactly the operation 
sequences counted by M. Let q,r he the number of consecutive I 1 O 2 or I 2 O 1 steps and 
number of times in the procedure given by s that the deque contains exactly one element 
respectively. So the contribution of s to M{a,u,x) is a^u'^x'^. 

Now we calculate the contribution to T from all the words of the form w, with s = 
X 1 X 2 ■ ■ ■ X 2 m fixed. Call this the contribution of s to T. Each Wi which begins (and ends) 
at height 1 in the bi-coloured Dyck path corresponding to w can be any tsip-word, whose 
steps from height 0 in its bi-coloured Dyck path are all red. Hence the possible words Wi 
are exactly those counted by Qi. Each Wi which begins and ends with the deque containing 
more than one element can be any tsip-word, and these are enumerated by Q. Also, the 
possible words v are exactly those counted by T. Since the half-length of rc is m plus the 
sum of the half-lengths of these sub-words, the contribution of s to r(l, n, 1) is 

u, 1YQ{1, u, 1). 

Now, each vertex at height 1 in the bi-coloured Dyck path corresponding to w occurs in 
one of the words Wi counted by Qi, or in v. The number of these vertices in any word 
Wi counted by Qi, is equal to the number of vertices in the bi-coloured Dyck path of Wi 
at height 0, which is 1 more than the power of x in the contribution of Wi to Qi{a,u,x). 
Since all other vertices at height 1 appear in v, the contribution of s to T{l,u,x) is equal 
to 

n™(xQi(l, u, x)YQ{l, u, x). 

Now we just need to consider the number of sub-words I 1 O 2 and I 2 O 1 in w. If we only 
consider the sub-sequences which occur in one of the words Wi or v, we would get the 
contribution 

u^{xQi{a, u, x)YQ{a, u, x). 

The only other situation where one of the these sub-sequences occurs is when one of the 
words Wi is empty, and the surrounding letters XjXj+i form one of these sub-words. Note 
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that this cannot happen at height 1, since w is top happy, so it only occurs in the case 
where Wi is counted by Q. This case is counted as 1 in Q{a,u) in the equation above, 
but it should be counted as a. Moreover, this is relevant to exactly q of the words Wi, 
exactly those where XjXj+i is either the word I1O2 or I2O1. Hence, the contribution of s 
to T{a, u, x) is 

vJ^{xQi{a, u, x)YQ{a, {Q(a, u) + a — l)'^T(a, u, x). 

Therefore, any top happy, unbreakable operation sequence s which contributes a'^u^x'^ to 
M{a,u,x), contributes 

to r(a, u, x), and this accounts for all of T(a, u, x) except for the 1 coming from the empty 
word. Hence, 

“ (1 + ^ f") ^ +1- 

□ 


Lemma 3 . 8 . The generating function D satisfies the equation 

D{t) = M(^l-^,tPYl^^ + l. ( 8 ) 

Proof. Using Proposition 2.5, we know that D{t) is the generating function for canonical 
operation sequences, counted by half-length. Using Lemma 3.6, we see that D{t) is the 
generating function for words w which decompose as w = xiwi... W2m-iX2mV, where 

• The (unbreakable) operation sequence s = X1X2 ■ ■ ■ X2m is top happy, 

• Each tsip word Wi is standard and outputs eagerly, 

• u is canonical, 

• If some XiXi+i is either I 1 O 2 or / 2 O 1 , then Wi is non-empty. 

As in the proof of the previous lemma, we will consider the contribution of any given 
top happy, unbreakable operation sequence s = X1X2 ■ ■ ■ X2m to D, assuming that s is 
counted by the monomial a^u^x^ in M(a,u, x). For each i, the word Wi can be any 
standard tsip word which outputs eagerly, except that it can’t be empty if XiXi+i is I1O2 
or I2O1 and there are q such values of i. Also, v can be any canonical operation sequence. 
Recall that standard tsip words which output eagerly are counted by P. Therefore, the 
contribution of s to D(t) is 

Hence, 

= ^ + 1 . 

□ 
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Theorem 3.9. Let D{t) he the length generating function for permutations which are 
sortable by a deque, and let P{t) be the length generating function for permutations which 
are sortable by two stacks in parallel. Then D and P satisfy the following two (equivalent) 
equations: 


P{t) 


2t{D -l-Dt) 


(9) 


and 

2D{f) = 2 + t + 2Pt- 2Pt^ - tVl - 4P + 4P2 _ gPH + 4.PH‘^ - APt. (10) 


Proof. First we substitute (6) into (7), to remove Qi: 


T{a, u,x) = M [ 1 + ^ ^ ^ ,uQ^, 


2x 


2Q — xQ + X J Q 


T 

+ 1 - 


Now we combine this with (4) to get 


T 

Therefore, 


1 


- 1 , 




p ’ (1 - 2py 


1 


, X = M 1 — —,tP, 


P' ' 2P-1-Px + x J 2P-1 


+ 1 . 


T - - 1 




P ’ (1 - 2P) 


2 ’ ■ 


P 


2 - - =M 1- -AP\l 


P’ 


2P- 1 


+ 1 . 


Using (8), we can write M in terms of D and P. So 


r(i-i o i\ (P-i)P T 

' P ’ (1 - 2P)2 ’ pj D 2P-1 


Solving for T gives the relation 


1 tP^ 


2 - 1 ^= ( 2 ^- 1 )^ 


P ’(1-2P)2’ Pj dP + P-D' 

Finally, using the formula (5) for T and rearranging gives the desired result. 


( 11 ) 

□ 


4 Analysis 

The main purpose of this section is to reduce the problem of showing that the generating 
functions P and D have the same radius of convergence to a few conjectures about the 
generating function Q{a,u) for quarter plane loops. We will also give some evidence for 
the stronger conjecture below. 

Conjecture 4.1. For n G Z>o, let pn be the number of permutations of size n which are 
sortable by two stacks in parallel and let dn be the number of permutations of size n which 
are sortable by a deque. We conjecture based on theorem 3.9 that pn ~ const ■ pp ■ and 
dn ~ const ■ pP ■ for some constants p. and 7. 

We first list three conjectures which are needed to show the above conjectures. The 
following are conjectures 10, 11 and 12 in [2], respectively. 
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Conjecture 4.2. The series Q{a,u) is (a + l)-positive. That is, Q takes the form 


Q{a,u) = '^u^Pnia + 1), 

n>0 


where each polynomial Pn has positive coefficients. 

Conjecture 4.3. the radius of convergence PQ{a) ofQ{a,-) is given by 


PQ{a) 


1 

, (2 + V2T^j2’ 

—a 


2(a-l)2’ 


if a >-112, 
if a e[-l,-1/2]. 


Conjecture 4.4. The series Qu{a,u) = ^ is convergent at u = pqia) for a > —1/3. 

Let tc denote the radius of convergence of the generating function P{t). Theorem 23 
in [2] states that assuming the two Conjectures 4.2 and 4.4 are both true, 


- 1 ), 

p(t)) 

for 0 < t < tc, with equality if and only if t = tc, and 

P{tc) < 


Moreover, the following corollary in [2] states that assuming the last three conjectures are 
all true, the following equation holds: 


\/2P{tc) — 1 + \/2tcP{tc). 


( 12 ) 


Note that in [2] the symbol S is used instead of P to denote the generating function for 
permutations sortable by two stacks in parallel. 

Theorem 4.5. Assuming the conjectures j.2, j.3 and j.j, the number of permutations 
of size n sortable by two stacks in parallel and the number of permutations sortable by a 
deque of size n have the same exponential growth rate. 

Proof. It suffices to prove that the generating functions P{t) and D{t) have the same radius 
of convergence tc. Since every permutation which is sortable by two stacks in parallel is 
also sortable by a deque, the coefficients of D{t) are no smaller than the coefficients of 
P{t). Hence the radius of convergence tjo of D(t) satisfies to < tc- Therefore, it suffices 
to prove that D{t) is convergent for t £ [0,tc)- 

Since P{t) is convergent for t G [0,tc)j the series 

1 - 4P(t) + 4P(t)2 - 8P{tft + 4.P{tft‘^ - 4:P{t)t 


is also convergent in this region. Since P{t) has positive coefficients, it is increasing on 
the interval [0,tc]- Hence, 

1 = H(0) < P{t) < Pitc) < 2 
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inside this interval. Therefore, since we are assuming Conjecture 4.3, we have the inequal¬ 
ity 

(2 - = (2 + ^2/P(i))2’ 

where equality holds if and only if t = tc- Taking square roots on both sides and rearrang¬ 
ing, noting that 2(2 — 1/P(t)) = [2 -\- y^2/P(t))(2 — y^2/P(t)), we get the inequality 

v/2P(t) > 1 + ^/2tP{t), 

with equality if and only if t = tc- Now we can remove the square roots as follows to get 
the inequality 

0 + 1 - V^){V^ - I + + l + V^) 

=1 - AP{t) + 4:P{tf - 8P{tft + APitft"^ - 4P{t)t, 

with equality if and only if t = tc- Since the series 

1 - 4P(t) -h 4P(t)2 - 8P{tft -h 4P(t)2t2 _ 4P(t)t 

is convergent and positive for t £ [0,tc)) the series 

V^l - 4P(t) + 4P(t)2 - 8P{t)H + 4P(t)2t2 _ 4P(i)t 

is also convergent in this domain. Hence, the series 

2D{t) = 2 + t + 2Pt- 2Pf - t\/l -4:P + AP"^ - 8PH + APH'^ - 4Pt 

is also convergent in this domain. □ 

For the following analysis, we will assume that 

P{t) = Co -h ci(l - t/tc) + Ca{l- t/tcT + o((l - t/tcT), (13) 

for some constants cq, ci and Cq, and a, with Cq ^ 0. In the next section we will present 
numerical evidence that a ~ 1.473, but for the following theorem we will assume only that 
2 > a > 1. 

Theorem 4.6. Assuming that P takes the form given in (13), and assuming the three 
conjectures j. 2 , j.3 and 4-4 j we have D{t) = D{tc) + A:d(1 — tftc)^^^ + o((l — t/tc)^/"^), 
where ko is given by the equation 

kD = -tc\l{“^coY/Hc - 2Cl^/^c. 

Proof. First, using (10), we can rewrite (12) as 

\/2co = 1-1- \/2tcCo- 


It follows that 

0 = (1 -|- y/2tcCo — \/2wo){l + y/2tcCo + \/2xf)){l — \/2tcCo — \/2w(j){l — \/2tfcQ + y/2wo) 

= 1 - 4co -F 4co - Scgtc -H 4cotc - 4cotc 
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Now, recall that 

2D{t) = 2 + t + 2Pt- 2Pt^ - 

Using (13) we can expand the expression under the square root as a power series in 
(1 - t/tc). 

1 - 4P + - SP'^t + 4P^t^ - APt 

=(1 - 4co + 4co - Scgtc + 4cotc - 4cotc) 

+(—4ci + ScqCi — Ihcocitc + Scocit^ — 4citc + Scptc + 4cotc — 8cQt^)(l — t/tc) 

+ (-4Cq, + ScoCq, - IhCoCatc + ScoCatl - ACatc){l “ t/tc)" 

+o((l - t/tc)") 

=(—4ci + ScqCi — Ihcocitc + Scocit^ — 4citc + Scotc + 4cotc — 8cotc)(l — t/tc) 

+ (-4Cq, + ScoCq, - IhCoCotc + ScoCat^ - 4cQ,tc)(l - t/tc)" 

+o((l - t/tc)") 

=gi(l - t/tc) + qa{l - t/tc)^^ + o((l - t/tc)"). 

Here qi and qa are constants defined by 

qi = —4ci + ScqCi — Ihcocitc + ScoCitc — 4citc + Scgtc + 4cotc — Scgtc 

and 

qa = —4Cq + ScqCq — IhCoCatc + ScoCatc — ACatc- 

Since this expression is non-negative for t G [0,tc], we must have qi > 0. 

Taking the square root of this, we get 

\/l - 4P + 4P2 _ 8PH + APH"^ - APt 

= (1 - t/tc)^^‘^Vqi + Qai'i- - t/tc)""^ + o((l - t/tc)"-l) 

= - t/tc)^/^ + 0((1 - t/tc)^/^) 

Finally, we can use this to determine the asymptotics of D 

2D{t) =2 + t + 2Pt- 2Pf - t Vl -AP + AP^ - 8PH + APH^ - APt 

=(2 + tc + 2cotc - 2cotc) - tcV^(l - t/tc)^/^ + o((l - t/tc)^/^) 

Simplifying this using the equation t/2co = 1 + \/2tcCo gives the desired expression for 
ko- □ 

Theorem 4.7. With the same assumptions as in the previous theorem, and the additional 
assumptions that ko > 0 and a < 3/2, we have the expansion 

D{t) = D{tc) + knil - t/tc)^/^ + ka{l - t/tc)"“^/^ + o((l - t/tc)"“^/^), 

where ka is given by the equation 


v/(2co)3/2tc - 2ciV^ 
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Proof. From the proof of the previous theorem we have 

\/l - 4P + 4P2 _ 8PH + 4 P 2 t 2 _ 

= (1 - tftc)^^‘^\/qi + Qai'i- - + 0((1 - 

Since qi > 0, we can expand this as a power series in (1 — t/tf) as follows 

(1 - t/tcY/‘^yJqi + qa{l - + o((l - 

=(i - + o((i - t/tcT-^)^ 

= ^(1 - + o((l - 

2 \/^ 

Therefore, we have 

2D{t) =2 + t + 2Pt- 2Pt^ - t\/l -AP + 4 P 2 _ _ 4 ^^ 

=(2 + tc + 2cotc - 2cotc) - + o((l - 

2V9i 

Simplifying this using the equation y/2co = 1 + ^/2tcCQ gives the desired expression for 
ka- □ 


5 Asymptotics 


In this section we make a numerical study of various generating functions related to 
quarter-plane loops, tsips and deques. We use the two most common methods of se¬ 
ries analysis, the ratio method and the method of differential approximants. Full details 
of these methods can be found in, for example, [ 6 ]. Both methods aim to estimate the 
radius of convergence (r.c.) Zc and associated exponent 6 of functions whose asymptotic 
behaviour is given by 


It follows that 


F{z) ~ A ^1 — —^ , as z ^ z^ . 

(14) 

A.r,-d-l 

(15) 


5.1 Ratio method. 

The ratio method, as the name implies, relies on extrapolating the ratio of successive 
coefficients, One has 

rn = -= — 1 -^ o(l/n) . (16) 

Jn-l Zc\n ) 

Clearly, extrapoltating the ratios against 1/n should, for sufficiently large n, give a 
linear plot which extrapolates to Ijzc at 1/n = 0. The gradient is — (0 -|- l)/zc. So from 
the intercept one can estimate the radius of convergence, and from the gradient and the 
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estimate of the radius of convergence, one can estimate the exponent 6. In favourable 
cases, where the singularity is precisely as given in (14), the term o(l/n) can be replaced 
by 0(l/n^), and the ratio plot will usually be linear from quite low values of n. 

If the r.c. is known, or very accurately estimated from some other method, it follows 
from (16) that a more precise estimate of the exponent can be made from estimators 

9n = n{l-Zc-rn)-l + o{l), (17) 


where in favourable cases the term o(l) can be replaced by 0(l/n). 

Even if the r.c. is not known, one can obtain an estimate of the exponent independent 
of the r.c. by extrapolating ratios of ratios, so that 


t 


n 


Tn 

1 



So we can define 0^^, an estimator of the exponent 0, as 


On = {tn - -1 = 0 + o(l). 


(18) 


(19) 


where again, in favourable cases, the term o(l) can be replaced by 0(l/n). 


5.2 Method of differential approximants 

The method of differential approximants [6] hts the known coefficients of a power series 
to a number (typically 10 or 12) of holonomic differential equations, and uses the critical 
parameters (the radius of convergence and exponent at that point) of those differential 
equations as estimators of the corresponding quantities for the underlying series expansion. 

More precisely, one uses the known series coefficients to find polynomials Qk{z) and 
P{z) such that the power series solution F(z) of the holonomic differential equation 



F{z) = F{z) 


( 20 ) 


agrees with the known coefficients of the function F{z) being approximated. The order 
M of the ode we refer to as the order of the approximant. 

Constructing such differential approximants (DAs) is straightforward computationally, 
and only involves solving a linear system of equations. Several such DAs are constructed 
by varying the degrees of the polynomials Qk{z) and F(z), while still using most, or all 
of the known series coefficients. The singularities are given by the zeros Zi, i = 1,, Nm 
of Qm{z), where Nm is the degree of Qm{z). We take as the dominant singularity that 
which is both closest to the origin and common to all (or almost all) the DAs. Critical 
exponents 6i follow from the indicial equation of the DA. For the simplest (and most 
frequent) situation where there is a single root of Qm{z) at Zi, 


ei = i-A/ + %hw' 

ZiQMiF) 

Slightly more complicated expressions are known for the cases of double, triple etc. roots. 
Further details of both methods can be found in [6] . 
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5.3 Quarter-plane loops 

We first studied quarter-plane loops, introduced by Albert and Bousquet-Melou in [ 2 ] and 
described in the previous sections. We generated 500 terms in the ogf from the recurrence 
relation given in the next sub-section. 

Using the method of differential approximants, in particular third-order approximants, 
we were able to estimate the critical point to an accuracy of between 6 and 20 significant 
digits, depending on the value of the parameter a. Based on these numerical results, we 
conjectured the a-dependence of the radius of convergence, which is given as Conjecture 
11 in [2], and is stated here immediately below Conjecture 4.1. For a = —1 and a = 1 
the ogf is holonomic, and we found the defining differential equation explicitly. That is 
to say, in those cases the differential approximants were precisely the defining odes. It is 
for these two values of a that 15 digit accuracy in the radius of convergence was obtained. 
For other values of a, the precision of our estimates was lower, which incidentally is good 
heuristic evidence that the underlying generating function is non-holonomic. 

We also confirmed the variation of the critical exponent with the parameter a, as 
reported in [ 2 ]. As a point of clarification, the exponents we refer to relate to the generating 
function, not the coefficients. For example, at a = 0, Uc = ( 6 -|-4\/2)“^, and we conjecture 

Q{u,0) = qo + qi{l - u/ucf + o((l - ujucf), 
whereas \u^]Q(u,0) ~ const ■ u~^ ■ where 9 = - 7 and qn and oi are 

^ ^ arccos(v2—1) 

constants. Similarly, we find 

Q (^U, -0 = go + gi(l - u/Ucf^'^ + o((l - u/Ucf/'^), 


where go and gi are constants, and Uc = 1/9, so that [u'^]Q{u, —1/2) ~ const ■ u~^ ■ n~'^l^. 

While we were unable to find an exact expression for the critical exponent for all values 
of a, we estimated the exponent values for several values of a > —1/2, and found that the 
exponent appears to be discontinuous at a = —1/2, which is also the case for half-plane 
and full-plane walks. In particular, we found that the exponent increased monotonically 
with a for a > —1/2. However it appears that the exponent when lim^_j._^/ 2 + is 1, but 
at a = —1/2 it is 3/4. Such monotonic, continuously varying exponents would preclude 
holonomic generating functions, except at isolated values of a. Subsequently, in as yet 
unpublished work, Kilian Raschel told us of his conjectured result for the exponent, which 
is 

-^ for a > 0. 


arccos 


a—1 


ct-f- 1 -f-U2T2I1 


Raschel’s conjecture agrees with our numerical results for an even broader range, notably 
a > —1/2. As mentioned, at o = —1/2 our series analysis gives 0.750 ± 0.002, from which 
we conjecture the exact value of 3/4. So not only do we observe the unusual phenomenon 
of a critical exponent varying steadily with a parameter, but also the phenomenon of a 
jump-discontinuity at a particular value of a, in this case at a = —1/2. 
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5.4 Deques and two stacks in parallel. 

Recall that D[t) is the ogf of the number of permutations sortable by a deque, and P{t) 
is the corresponding ogf for the number of permutations sortable by two stacks in parallel 
(tsips). Then, as shown above in Theorem 3.9, 


^ ’ 2t{D{t) - 1 - tD{t)) ’ 

and equivalently that 

D{t) = ^-+l+tP{t)-t^P{t)-'^-^/l - AP{t) + 4P(t)2 - 8tP(t)2 + 4 t 2 p(t )2 
The series start: 


( 21 ) 


4tP(t). (22) 


P[t) = l + t + 2f + Qt^ + 23t^ + 103t^ + 513t® + 2760t^ + 15741t® + • • • tsips, 
D{t) = l + t + 2t^ + + 24:t^ + + 634t® + 3762t^ + 23638t® + • • • deques. 


We generated 500 terms in the tsip series from the functional equations in [2] using the 
following method: 

For integers n, fc, x, y, let s(n, k, x, y) denote the number of n step quarter plane walks 
with k ES or NW corners which start at (0, 0) and end at (x, y), so that 

OO OO 

Q{a, u) = s(n, k, 0, 0)a^u”. 

71=0 k=0 

We calculated the values s(n, k, x, y) for n < 500 via the recurrence relation 

s(n, k, X, y) =s{n — 1, k, x — l,y) + s{n — l,k,x,y — 1) 

+s(n - 1, fe, X + l,y) + s{n - l,k,x,y + 1) 

+s(n — 2,A: — l,x + l,y — 1) — s{n — 2,k,x + l,y — 1) 

+s{n — 2, k — l,x — l,y + l) — s(n — 2, A:, x — 1, y + 1). 


Now that we have calculated coefficients of Q{a,u), we can expand 


F{p,t) = Q (^-p, 


2y+ 1 


as a series in p and t. Next we calculated the first 500 terms of the series p{t) which sends 
F{p{t),t) to 0. Note that p is called S* in [2]. Finally we calculated the first 500 terms of 
P{t), using the equation 


P{t) 


1 

1-pit)' 


Calculating the series P{t) from the terms in Qia, u) is very fast, so most of the time in 
this algorithm was spent calculating the coefficients of Q. 

Using the tsip series, along with Theorem 3.9, we then generated the first 500 terms 
in the deque series. We subjected these two long series to differential approximant [ 6 ] 
analysis. 
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For tsips, 6th order DAs show the dominant singularity to be at tc = 0.1207524975763(2) 
with an exponent 1.47309(3), with another exponent with the value 1.94652(2) at tc = 
0.1207524975763(1). A third exponent at the same place with the value 4.72(4) is also 
suggested. The quoted errors reflect only the scatter in individual approximant estimates, 
and because of the presence of confluent singularities should be multiplied by a factor of 
10 at least to be on the safe side. 

With 8th order DAs, we And the dominant singularity to be at tc = 0.1207524975763(2) 
with an exponent 1.47309(4), and with another exponent 1.94652(2) at tc = 0.1207524975764(8). 
A third exponent at the same place with the value 4.72(4) is also suggested. 

For the deque series, the unbiased analysis showed the dominant singularity to be at 
Zc = 0.12075249773(4), with exponents values of 0.507(4), 0.970(2), and 1.414(1). The 
uncertainties quoted just reflect the variability in the estimates across many DAs, and 
the lack of overlap between Zc estimates for deques and tsips beyond the 10th significant 
digit suggests that they are too optimistic. Nevertheless, 10-digit agreement gives consid¬ 
erable credence to the conjecture that they are indeed equal, and we assume this for the 
subsequent analysis. The dominant exponent is very close to 1/2 exactly, which provides 
numerical support for the conjectured square-root singularity obtained in the previous 
section. 

Another way to study the exponents for these two problems is by the ratio method. 

We write [t”]T’(t) ~ const ■ and [t^]D{t) ~ const ■ . Then from the Hadamard 

(coefficient-by-coefficient) quotient, 

e 

Qn = 7 — , ^ ' ~ const ■ n , 

^ [t^]D{t) 

where 6 = gp — Qd- Then the ratios of successive coefficients rn behave as 

r„ = ~ 1 -I- - -h o(l/n). 

qn-i n 

We can estimate the exponent 9 by extrapolating a sequence of estimators {On}, defined by 
n ■ {cn — 1) ~ + o(l)- The result is shown in Figure 4, in which 9n is shown extrapolated 

against 1/n. While it is difficult to extrapolate this curve, a limit in the range [-0.975- 
-0.972] looks plausible. 

We also analysed the deque series by the ratio method. From eqn (19), estimators of 
the deque exponent can be found, and these are shown extrapolated against 1/n in Fig 
5. The plot is quite linear and is clearly going to a value close to —1.5. 

We can refine this by calculating the exponent g more precisely from the local gradient 
of the previous curve. This is hn = {gn — gn-i)n{'y — n), where gn is the estimator of 
gd- Then a straight line with this gradient will meet the ordinate at gn — which should 
be a refined estimator of the exponent g. This plot is shown in Fig 6. That the curve 
appears to be going slightly below -1.5 is, we believe, of no consequence. We believe that 
if we had several thousand terms we would see this curve pass through a minimum, and 
increase to -1.5 exactly. 

Accepting the conjectured square-root singularity of the deque generating function, 
numerical results thus far suggest that 

D{t) = D{tc) + Di{tc)yjl-t/tc + D2{tc) ■ (1 - t/tcf/^^^ + ■■■ (23) 


25 





0 0.005 0.010 0.015 0.020 



Figure 4: Exponent estimates On plotted against 1/n. Note limit appears to be around 
-0.974. 



2 

n 


Figure 5: Exponent estimates Qn plotted against 1/n for deques. The limit appears to be 
5Ri-1.5. 
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Figure 6: Extrapolated exponent estimates plotted against 1/n for deques. The limit 
appears to be around g ^ —1.500. 

where 0 < A ~ 0.47 < 1/2, as this gives the value 0.97 observed as the sub-dominant 
exponent for deques. Similarly, it appears that 

P{t) = P{Q - tcP'{tc){l - t/p) + Pi{tc){l - + o((l - t/tcf+^). (24) 

To see how this is consistent with (21) and (22), write (23) as 

D{t) = D{tc) + Di{tc)il - + D2{tc){l - t/tcf H-, 

where /? = 1/2-|-A. Substitute into (21). This gives an expression for P{t) which includes 
terms of 0{y/l — t/tc) and 0((1 — Then a little algebra shows that the coefficients 

of both these terms vanish if 

1 + 

= Tviy (25) 

Remarkably, this is the case, as follows from the results of (12). 

More precisely, in the proof of (4.6) we have shown that 

= 2(l-Vt/)2' 

and 

= 2 + tc P 2P{tc){tc — t^), 

which can readily be shown to give (25). 

From eqn (24) and the proof of (4.6), we have 

Di{tc) = -2^/4 • + ^/P{tc) • tc{l - V^c) • P^{tc). (26) 
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We also estimated the value of P' (tc) numerically from Fade approximants to the series 
for P{t) — P{tc), evaluated at t = tc- Numerical estimates of P'{tc) give Di{tc) ~ —0.0540 
from (26), which is precisely equal to direct estimates of di(fc) obtained from our numerical 
analysis of the deque series. There we formed the series for 

D{t) - D{tc) 

- t/tc 

and evaluated the approximants att = tc- Furthermore, using this amplitude value Di{tc), 
and subtracting D{tc) + Di{tc)\/l — t/tc from the deque series D{t) gives a remainder 
series that behaves as const ■ (1 — t/tcY, where 6 ~ 0.973. So both ratio and differential 
approximant analyses clearly identify this confluent exponent. 

Similarly, we analysed the series for two stacks in parallel by the ratio method. The 
exponent estimators gn, given by eqn (17), are plotted against 1/n in Fig. 7. The plot is 
visually quite straight and is clearly going to a limit of about -2.475, as shown in Fig 7. 



]_ 

n 


Figure 7: Exponent estimates Qn plotted against 1/n for two stacks in parallel. The limit 
appears to be about g ^ —2.475. 

As with the deque series, we can refine this exponent estimate by calculating the 
exponent g more precisely from the local gradient of the previous curve. The refined 
estimate of the exponent g is shown in Fig 8, and from that plot a limit around -2.474 
seems quite plausible. 

So the ratio plots support the conclusion that the exponent for the deque generating 
function is very close to 0.5, corresponding to a square-root singularity, as conditionally 
proved in Theorem 4.6. The first confluent exponent for deques has the approximate value 
0.973, and these two exponents add to give the observed value 1.473 for the exponent of 
the ogf for tsips, as can be seen from (21). A direct ratio analysis of the ogf for tsips also 
gives an exponent value around 1.472 — 1.474. 
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Figure 8: Extrapolated exponent estimates gn plotted against 1 /n for two stacks in parallel. 
The limit appears to be around g ~ —2.474. 


Finally, we calculate the amplitudes, or pre-multiplying constants Kd and Kp, where 

[t'^]D{t) = dn ^ Kd ■ and [t'^]P{z) = Pn ^ i^p ■ 

We do this by forming the sequences dn-tc ■ ~ Kd + o(l) and Pn'tc ' ~ /tp + o(l), 

and extrapolating these against 1/n. In this way we estimate Kd = 0.01524 ± 0.0005, and 
Kp = 0.08025 ± 0.0010. 

Moving now from the arena of careful numerical work to that of wild speculation, in 
some unpublished work we have studied the behaviour of 421-3 pattern-avoiding permu¬ 
tations. We identihed the critical exponent numerically, as 

^ (^1 = 1.472799717437.... 

We were struck by the fact that this is tantalisingly close to the observed exponent of 
two stacks in parallel, though there is no obvious reason that the two problems should be 
connected, except that they both involve pattern-avoiding permutations. 

6 Conclusion 

For quarter-plane loops we have provided numerical support for the variation of the critical 
point with corner-parameter a, as conjectured in [2]. We also estimated the value of the 
critical exponent for a variety of values of a > —1/2. Our numerical exponent values 
agree with the conjectured formula provided, in private correspondence, by Kilian Raschel. 
There appears to be an exponent discontinuity at a = —1/2. 
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Our principal result is Theorem 3.9, giving the solution of the deque generating func¬ 
tion in terms of that for tsips. Other conclusions are subject to the validity of certain 
conjectures. These include the result that the critical point for both the deque and tsip 
generating functions are equal. We also provide compelling numerical evidence for this, 
finding estimates that agree to 10 significant digits. Subject to these conjectures, we prove 
in Theorem 4.6 that the deque generating function has a square-root singularity. 

From the solutions of the deque and tsip generating functions, we produced 500 terms 
of the generating functions, and subjected these to careful numerical analysis. The asymp¬ 
totic form of the generating functions was found to be, for deques, 

D{t) ~ D{tc) + Di{Q^l-t/tc + D2{tc) • (1 - + • • • (27) 

3/2 

where A « 0.473, 4 = 0.1207524977, do(4) = « 1.185059767, and diiQ ^ 

—0.0543. For tsips we found 

Pit) ~ + -tcP' (tc) (1 - t/tc) +Plitc)il- t/tc) + P2 (tc) (1 - t/tc) + • • • , (28) 

where T’(tc) = ~ 1.174361446, T’i(tc) ~ 0.1940, and A ss 0.473. At the coefficient 

level we have 

[t^]Diz) « 0.01524 • t-^ ■ 

and 

[r]P(t) PS 0.08025 • 

Given the simplicity of the relationship (1) between the two generating functions Pit) 
and Dit) characterising tsips and deques respectively, there could be a much simpler proof 
than the one we have constructed. It also remains a source of some frustration that we 
cannot establish the asymptotics rigorously, nor even unequivocally prove that the critical 
points of deques and tsips are identical, but hopefully our results will stimulate work in 
this direction. 
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